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Quantum gravity in an AdS spacetime is described by an SU{N) Yang-Mills theory 
on a sphere, a bounded many-body system. We argue that in the high temperature phase 
the theory is intrinsically non-perturbative in the large limit. At any nonzero value 
of the 't Hooft coupling A, an exponentially large (in A^) number of free theory states 
of wide energy range (or order A) mix under the interaction. As a result the planar 
perturbation theory breaks down. We argue that an arrow of time emerges and the dual 
string configuration should be interpreted as a stringy black hole. 
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1. Introduction 

While tlie equations of general relativity are time symmetric themselves, one often 
finds solutions with an intrinsic arrow of time, due to the presence of spacelike singular- 
ities. Familiar examples include FRW cosmologies and the formation of a black hole in 
a gravitational collapse. In the case of a gravitational collapse to form a black hole, the 
direction of time appears to be thermodynamic, since a black hole behaves like a thermo- 
dynamical system [|T],|^,0] . It has also long been speculated that the thermodynamic arrow 
of time observed in nature may be related to the big bang singularity . 

In an anti-de Sitter (AdS) spacetime, a microscopic understanding of the emergence of 
thermodynamic behavior in a gravitational collapse can be achieved using the AdS/CFT 
correspondence which states that quantum gravity in an asymptotic AdS^ x 

spacetime is described by an A/" = 4 SU{N) super- Yang-Mills (SYM) theory on an 

The classical gravity limit of the AdS string theory corresponds to the large N and 
large 't Hooft coupling limit of the Yang-Mills theory. A matter distribution of classical 
mass M in AdS, can be identifiedB with an excited state of energy E = ^N"^ in the SYM 
theory with /i a constant independent of A^. The gravitational collapse of the matter 
distribution can be identified with the thermalization of the corresponding state in SYM 
theory, with the resulting black holei identified with thermal equilibrium P,^. In this 
context, it is natural to suspect that the appearance of a spacelike singularity at the end 
point of a collapse should be related to certain aspect of thermalization in the SYM theoryl 

A crucial element in the above description is the large limit. Af = 4 SYM theory on 
is a closed, bounded quantum mechanical system with a discrete energy spectrum. At 
any finite A^, no matter how large, such a theory is quasi-periodic (i.e. has recurrences), time 
reversible, and never really thermalizes. However, to match the picture of a gravitational 
collapse in classical gravity, an arrow of time should emerge in the large A^ limit for the 
SYM theory in a generic state of energy E = fj.N'^ with a sufficiently large /i. This 
consistency requirement immediately raises several questions: 
1. What is the underlying physical mechanism for the emergence of an arrow of time in 

Yang-Mills theory? 

^ See Appendix A for a brief review of parameter translations in AdS/CFT. 

^ Assume M is sufficiently big that a big black holes in AdS is formed, which also implies that 
fi should be sufficiently big. 

^ Some interesting ideas regarding spacelike singularities and thermalization have also been 
discussed recently in []To| . 
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2. Is large 't Hooft coupling needed? 

3. Suppose an arrow of time also emerges at small 't Hooft coupling, what would be 
the bulk string theory interpretation of the SYM theory in such an excited state? A 
stringy black hole? Does such a stringy black hole have a singularity? 

4. Is there a large N phase transition as one decreases the 't Hooft coupling from infinity 
to zero? 

It would be very desirable to have a clear physical understanding of the above questions, 
which could shed light on how spacelike singularities appear in the classical limit of a 
quantum gravity and thus lead to an understanding of their resolution in a quantum 
theory. 

The emergence of an arrow of time is also closely related to the information loss 
paradoxal. At finite A^, the theory is unitary and there is no information loss. But in 
the large N limit, an arrow of time emerges and the information is lost, since one cannot 
recover the initial state from the final thermal equilibrium. Thus the information loss 
in a gravitational collapse is clearly a consequence of the classical approximation (large 
N limit), but not a property of the full quantum theory. While AdS/CFT in principle 
resolves the information loss paradox, it remains a puzzle whether one can recover the 
lost information using a semi-classical reasoning. From this perspective it would also be 
valuable to understand the various questions listed in the last paragraph. 

The purpose of the paper is to suggest a simple mechanism for the emergence of an 
arrow of time in the large N limit and to initiate a statistical approach to understanding 
the quantum dynamics of a YM theory in highly excited states. In particular, we argue 
that the perturbative planar expansions break down for real-time correlation functions and 
that there is a large N "phase transition" at zero 't Hooft couplingi. We also argue that 
time irreversibility occurs for any nonzero value of the 't Hooft coupling. 

The plan of the paper is as follows. In section 2 we introduce the subject of our study: 
a family of matrix quantum mechanical systems including A/" = 4 SYM on . We highlight 
some relevant features of the energy spectrum of these theories. Motivated by the classical 
mixing properties, we introduce observables which could signal time irreversibility. The 
simplest of them are real-time correlation functions at finite temperature, which describe 



^ To our knowledge this connection was first pointed out in the context of AdS/CFT in 

^ See e.g. [|, 11, 12, 13, 14, 15. If;, 17| for recent discussions. 
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See [18, 19 1 for some earlier discussion of a possible large A'' phase transition in A. 
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non-equilibrium linear responses of the systems. The rest of the paper is devoted to 
studying these observables, first in perturbation theory, and then using a non-perturbative 
statistical method. In sec 3 we compute real-time correlation functions in perturbation 
theory. We find that at any finite order in perturbation theory, the arrow of time does 
not emerge. In sec 4 we argue that the planar perturbative expansion has a zero radius of 
convergence and cannot be used to understand the long time behavior of the system. In 
section 5 we give a simple physical explanation for the breakdown of perturbation theory. 
We argue that for any nonzero 't Hooft coupling, an exponentially large (in A^^) number 
of free theory states of wide energy range (or order A^) mix under the interaction. As a 
consequence small A and long time limits do not commute at infinite A^. In section 6 we 
develop a statistical approach to studying the dynamics of the theories in highly excited 
states, which indicates that time irreversibility occur for any nonzero 't Hooft coupling A. 
We conclude in section 7 with a discussion of implications of our results. 

2. Prelude: theories and observables of interest 

In this section we introduce the systems and observables we want to study. 

2.1. Matrix mechanical systems 

We consider generic matrix quantum mechanical systems of the form 

S = mrjdt ^ Q(AMJ2 _ ^ulM^^ - j dt V{M^- A) (2.1) 

which satisfy the following requirements: 

1. Mq, are N x N matrices and DtM^t = df — i[A, M^] is a covariant derivative. One can 
also include fermionic matrices, but they will not play an important role in this paper 
and for simplicity of notations we suppress them. 

2. The frequencies uJa in ( p.l|) are nonzero for any a, i.e. the theory has a mass gap and 
a unique vacuum. 

3. The number of matrices is greater than one and can be infinite. When there is an infi- 
nite number of matrices, we require the theory to be obtainable from a renormalizable 
field theory on a compact space. 

4. V(Mq;A) can be written as a sum of single-trace operators and is controlled by a 
coupling constant A, which remains fixed in the large limit. 
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A/" = 4 SYM on is an example of such systems with an infinite number of matrices 
(including fermions) when the Yang-Mills and matter fields are expanded in terms of 
spherical harmonics on (see e.g. [pO| , |2T[] ). In this case, Wa are integer or half-integer 
multiples of a fundamental frequency uq = 1/R with R the radius of the S'^. The number of 
modes with frequencies u>a — ^ increases with k as a power. V(Mo,; A) can be schematically 
written as0 

V = N (yXVsiM^) + XV^{M^)^ (2.2) 

where V3 and V4 contain infinite sums of single-trace operators which are cubic and quartic 
in Mot and dtMa. A = Qym-^ the 't Hooft coupling. 

In this paper we work in the large limit throughout. Our discussion will only depend 
on the large scaling of various physical quantities and not on the specific structure of 
the theories in ( |2.1| ) like the precise field contents and exact forms of interactions. For 
purpose of illustration, we will often use as a specific example the following simple system 

S = — tr / dt [(AMi)2 + {DtM2f - ooliMl -f M|) - AM1M2M1M2] . (2.3) 

2.2. Energy spectrum 



( pTl| ) has a U{N) gauge symmetry and physical states are singlets of U{N). One can 
classify energy eigenstates of a theory by how their energies scale with in the large A^ 
limit. We will call the sector of states whose energies (as measured from the vacuum) 
are of order 0(1) the low energy sector. As motivated in the introduction, we are mainly 
interested in the sector of states whose energies are of order fj.N'^ with fj, independent of 
A^, which will be called the high energy sector. The density of states in the low energy 
sector is of order 0(1), i.e. independent of A^, while that of the high energy sector can be 
written in a form 

Q(E) ~ e''^^^^\ E = fiN^ (2.4) 



with s{n) some function independent of A^. (|2.4|) follows from the fact that the number of 
ways to construct a state of energy of order 0(A^^) from 0{N'^) oscillators of frequency of 
0(1) is an exponential in A^^. The presence of interaction should not change this behavior 



^ The precise form of the interactions depends on the choice of gauge. It is convenient to 
choose Coulomb gauge V • ^ = 0, in which the longitudinal component of the gauge field is set 
to zero. In this gauge, Ma include also non-propagating modes coming from harmonic modes of 
ghosts and the zero component of the gauge field. 
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at least for n sufficiently large. (|2.4| ) is the reason why we restrict to more than one matrix 
in ( pj.ll ). For a gauged matrix quantum mechanics with a single matrix one can reduce 
the matrix to its eigenvalues and ( p.4| ) does not apply. When /i is sufficiently large, 
should be a monotonically increasing functioni of // and we will restrict our definition of 
high energy sector to such energies. 

For A/" = 4 SYM, states in the low energy sector correspond to fundamental string 
states in the AdS spacetime, while the states in the high energy sectors may be considered 
as black hole microstateJl. 

A convenient way to study a system in excited states is to put it in a canonical 
ensemble with a temperature T = ^. The partition function and free energy are defined 
by (Tr denotes sum over all physical states and H is the Hamiltonian) 

Z = Tre-^^ = e-^^ . (2.5) 

We will always keep T fixed in the large limit. Below low and high temperature refers to 
how the temperature is compared with the mass gap of a theoryS. As one varies T, differ- 
ent parts of the energy spectrum are probed. For the family of matrix quantum mechanical 
systems (|2.1[) , there are two distinct temperature regimes. At low temperature, one probes 
the low energy sector and the free energy F is of order At high temperature F is of 

order 0(A^^) and the high energy sector is probed. It may seem surprising at first sight 
that one can probe the sector of energies of 0{N'^) using a temperature of 0(1). This is 
due to the large entropy factor ( |2.4| ) which compensates the Boltzmann suppression. For 
A/" = 4 SYM theory at strong coupling, there is a first order phase transition separating 
the two regimes at a temperature of order where R is the AdS radius p^,|6|,|8[| . A first 
order phase transition has also been found for various theories in the family of at 



weak coupling p3|j24,21 



An important feature of the high energy sector is that the large N limit is like a 
thermodynamic limit with A^^ playing the analogous role of the volume factor. In this 
limit the number of degrees of freedom goes to infinity while the average excitations per 
degree of freedom remain finite. The thermal partition function 

Z{(3) = Tre"^^ = j dEn{E)e-^^ (2.6) 



That is, the theory should have a positive specific heat for sufficiently large. 
^ Note that at a sufiiciently high energy, the most entropic object in AdS is a big black hole. 
For example for J\f = 4 SYM on S^, low (high) temperature means T ^ {T ^) 
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is sharply peaked at an energy ~ 0{N'^) (with a width of order 0{N)) determined by 

dS{E) 



= /3, S{E)=\ogn{E) (2.7) 



dE 

Note that the leading dependence of S{E) has the form S{E) = N'^s^y) (see (|2.4|)) with 
11 = E/N'^ characterizing the average excitations per oscillator degree of freedom. Equation 
( p.7| ) can also be interpreted as the equivalence between canonical and micro canonical 
ensemble0. Note that since F ~ 0{N'^), the high temperature phase can be considered a 



"deconfined" phase [25 



2.3. Observables 

In a classical Hamiltonian system, time irreversibility is closely related with the mixing 
property of the system, which can be stated as follows. Consider time correlation functions 

CAB{t) = {A{^'X)B{X)) - {A) {B) (2.8) 

where A, B are functions on the classical phase space parameterized by X. describes 
the Hamiltonian flow, where is a one-parameter group of volume-preserving transfor- 
mations of the phase space onto itself. (...) in (|2.8[) denotes phase space average over a 



constant energy surface. The system is mixingirl iff [26 



CAB{t)^Q, t^oo (2.9) 

for any smooth functions A and B. 

The closest analogue of ( |2.8| ) for the matrix quantum mechanical systems we are 
considering would be 

G,it) = (z|O(t)O(0)|z) - {i\O{0)\if (2.10) 



In contrast in the low energy sector, since both the free energy and the density of states are 
independent of A'^, in generic models there is no large parameter that one can use to perform the 
saddle point approximation to equate two ensembles. 

Note that mixing is a stronger property than ergodic which involves long time average. The 
ergodic and mixing properties can also be characterized in terms of the spectrum of the Koopman 
operator. For example, a system is mixing iff the eigenvalue 1 is simply degenerate and is the 
only proper eigenvalue of the Koopman operator [ p6| . 
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where \i) is a generic energy eigenstate in the high energy sector, and O is an arbitrary 
gauge invariant operator which when acting on the vacuum creates a state of finite energy 
of order More exphcitly, denoting li^o) = C'(O) with \Q) the vacuum, we require 

{ipc>\H\ipo) ~ Note that for A/" = 4 SYM on S^, a local operator 0{t,x) of dimension 

0(1) on 5""^ is not aUowed by this criterion since x) creates a state of infinite energy. To 
construct a state of finite energy one can smear the local operator over a spatial volume, 
e.g. by considering operators with definite angular momentum on S^. Without loss of 
generality, we can take O to be of the form 

O = tr(M,, ■ ■ ■ M,„^ )tr(M^, ■ ■ ■ J ■ ■ ■ tr(M^, ■ ■ ■ M^„^ ) (2.11) 

with the total number of matrices K = XliLi ^fc independent of A^. We will call such 
operators small operators. The reason for restricting to small operators is that they have 



a well defined large limit in the sense defined in [^. More explicitly, if one treats 
the large limit of a matrix quantum mechanics as a classical system, then (|2.11| ) with 
K ~ 0{1) are smooth functions on the corresponding classical pase space. From AdS point 
of view, such operators correspond to fundamental string probes which do not deform the 
background geometry. If for all small operators O and generic states |i) in the high energy 
sector 

Gi{t)^0, t^oo (2.12) 

one can say the system develops an arrow of time. In particular, (|2.12| ) implies that one 
cannot distinguish different initial states from their long time behavior (i.e. information is 
lost). 

Energy eigenstates are hard to work with. It is convenient to consider microcanon- 
ical or canonical averages of ( |2.10|) , for example, the thermal connected Wightman func- 
tions (see Appendix B.l for a precise definition of "connected" and the constant C below) 

G+{t) = m)Om^ = ^Tr (e-^^0(t)0(O)) - C (2.13) 

or retarded functions 

Gnit) = ^Tr{e-^^[0{t),Om) ■ (2-14) 

We shall take the temperature T to be sufficiently high so that Ep determined from ( |2.7| ) 
lies the high energy sector. Equation ( p.l2|) implies thatlll 

G'fl(t)^0, G+{t)^0, t^+oo. (2.15) 
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( 2.12 ) in fact implies the following to be true for any ensemble of states. 
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Note that Gfi{t) measures the hnear response of the system to external perturbations 
caused by O. That Guit) for t ^ oo imphes that any small perturbation of the 
system away from the thermal equilibrium eventually dies away. In a weaker sense than 
( p.l2|) , (|2.15| ) can also be considered as an indication of the emergence of an arrow of time. 



In frequency space, the Fourier transform^ of ( |2.13| ) and ( f2.14| ) can be written in 
terms of a spectral density function p{uj) (see Appendix B.l for a review) 



Gr{u) = 



d.' p(.') (2.16) 



OO 



2n u — u' + ie 



( p.l5|) may be characterized by properties of the spectral density p{uj). For example from 
the Riemann-Lebesgue theorem, ( |2.15| ) should hold if p{u)) is an integrable function on the 



real axis. Since other real-time correlation functions can be obtained from (or spectral 
density function p{uj)) from standard relations, for the rest of the paper, we will focus on 
G-i- only. 

For A/" = 4 SYM at strong coupling, it is convenient to take O to have a definite 
angular momentum I on S^. (|2.13|) and ( |2.14|) can be studied by considering a bulk field 



propagating in an eternal AdS black hole geometry and one does find the behavior ( |2.15| ) 
as first emphasized in 0. In the bulk language, ( p.l5|) can be heuristically interpreted 



as the fact that any small perturbation of the black hole geometry eventually dies away 
by falling into the horizon. Furthermore, by going to frequency space, one finds that the 
Fourier transform G_|_(t<;,/) has a rich analytic structure in the complex w-planelll, which 
encodes that the bulk black hole geometry contains a horizon and singularities. The main 
features can be summarized as follows |p8|| : 

1. G+{u!, I) has a continuous spectrum with u G (— oo, +oo). This is due to the presence 
of the horizon in the bulk. 

2. In the complex w-plane, the only singularities of G+{ujJ) are poles. The decay rate 
for G+ (t) at large t is controlled by the imaginary part of the poles closest to the real 
axis, which is of order (3. 



We use the same letter to denote the Fourier transform of a function, distinguishing them 
by the argument of the function. 

Similar things can also said about Gr{uj, I) which can be obtained from G+(a;, I) using stan- 
dard relations. 
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3. The presence of black hole singularities in the bulk geometry is encoded in the behavior 
of G_|_(a;, I) at the imaginary infinity of the c<;-plane0. In particular, 
3a. G+{u, I) decays exponentially as w — ±zoo. 

3b. Derivatives of (^+(0;, /) over / evaluated at / = are divergent as u; — > ±zoo. 
As emphasized in none of the above features survives at finite A^, in which case0 



has a discrete spectrum and is a sum of delta functions supported on the real axis. This 
indicates that concepts like horizon and singularities only have an approximate meaning 
in a semi-classical limit (large N limit). 

To understand the information loss paradox and the resolution of black hole singular- 
ities, we need to understand how and why they arise in the classical limit of a quantum 
gravity. In Yang-Mills theory, this boils down to understanding what physics is missed in 
the large N limit and why missing it is responsible for the appearance of singularities and 
the loss of information. With these motivations in mind, in this paper we are interested 
in understanding the following questions 

1. Can one find a qualitative argument for the emergence of an arrow of time in the large 
N hmit? 

2. Does the analytic behavior observed at strong coupling persist to weak coupling? 
which we turn to in the following sections. 



3. Non-thermalization in perturbation theory 



In this section we consider (|2.13|) in perturbation theory in the planar limit. We will 
find that real-time correlation functions have a discrete spectrum and oscillatory behavior. 
Thus the theory does not thermalize in the large N limit. 

In perturbation theory, Gj^{t) can be computed using two methods. In the first 
method, one computes Ge{t) with < r < /3 in Euclidean space using standard Feynman 



See also [ ]29| for signature of the black hole singularities in coordinate space. 

Note that even though TV = 4 SYM on is a field theory, at finite N the theory can be 
efi^ectively considered as a theory with a finite number of degrees of freedom, since for any given 
energy E, there are only a finite number of modes below that energy. Furthermore, given that the 
number of modes with frequency grows with k only as a power, it is more entropically favorable 
to excite modes with low k for E ~ O(Y^) and modes with cOa ~ 0{N) are almost never excited. 
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diagram techniques. (t) can then be obtained by taking t = it + e. An alternative way 
is to double the fields and use the analogue of the Schwinger-Keldysh contour to compute 



the Feynman function Gpi^^) in frequency space ||3^, from which G+{uj) can be obtained. 
In the Euclidean-time method it is more convenient to do the computation in coordinate 
space since one does not have to sum over discrete frequencies, while in the real-time 
method frequency space is more convenient to use. 
We look at the free theory first. 

3.1. Free theory 

To evaluate ( p.l3|) in free theory, it is convenient to use the Euclidean method. The 
Euclidean correlator 

Gg)(r) = (0(r)0(0))o,^, < r < /3 (3.1) 
with O of the form (|2TTD can be computed using the Wick contractionlll 



M-(r)M-(0) = ^ E 9'^E\r-m(3;u^jUru;- (3.2) 



^ ' m= — oo 



where g^^^ is the propagator at zero temperature 



In ( |3.2D U is a unitary matrix which arises due to covariant derivatives in (^TT|) and can be 
understood as the Wilson line of A wound around the r direction. In the evaluation of free 
theory correlation functions (■ ■ ■)q ^ in ( p.l| ), one first preforms the Wick contractions ( |3.2| ) 
and then performs the unitary matrix integral over U, which plays the role of projecting 
the intermediate states to the singlet sector. In the large N limit, the U integral can be 
evaluated by a saddle point approximation. Note in particular that 

[/ -> 1, T ^ oo (3.4) 

Equation ( p.4| ) indicates that the singlet condition should not play an important role for 
states of sufficiently high energy. 



18 



see e.g. |31] for a derivation of the following equation and some examples of correlation 



functions in free theory. 
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For definiteness, we now restrict to theories with a single fundamental frequency ujq 
like A/" = 4 SYM or ( |2.3| ). Relaxing this restriction does not affect our main conclusions, 
as will be commented on in various places below. Wick contractions in give rise to 
terms of the form e'^^°'^ for some integer n, while the [/-integral computes the coefficients 
of these terms. Thus ( |3.1| ) always has the form 



G^E\r)= E ^niP)^-'^- (3-5) 

Tl=-A 



where A is the dimension of the operatoiil^. Analytically continuing (|3.5| ) to real time, we 
find that 

A 

Gfit)= c^(/3)e-^"'^°* (3.6) 

n=-A 

and 

A 

G^iu;) =2n J2 Cn(/3)5(w - nuo) • (3.7) 

n=-A 

Thus in the large N limit, the correlation function always shows a discrete spectrum 
is quasi-periodic. The results are generic. If the theory under consideration has several 
incommensurate fundamental frequencies, one simply includes a sum like those ( p.5|) and 
( p.7|) for each such frequency. The maximal number of independent exponentials is 2^, 
where K is the total number of matrices in O. This is due to that each matrix in O can 
only connect states with a definite energy difference. 

It is also instructive to obtain (|3.6| ) using a different method. By inserting a complete 
set of free theory energy eigenstates in (|2.13|) we find that 



^0 , 

a,o 

where \a) is a free theory state with energy and pab = \ {O'\C^{0)\b) p. To understand 
the structure of ( |3.8| ) we expand O{0) in terms of creation and annihilation operators 
associated with each (Mq,)^^, from which we find that 
A. Due to energy conservation, O can connect levels whose energy differences lie between 
— Awo and Awq, i-e. pab can only be non- vanishing for |ea — etl < Auq. 



Note that for A/" = 4 SYM the dimension of Ma is given by For other matrix quantum 
mechanical systems without conformal symmetry one can use a similar definition in free theory. 
For bosonic operators, A are integers. 
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B. O can only connect states whose energy differences are integer multiples of ojq i.e. pab 
can only be non- vanishing for — = nujQ with |n| < A integers (or half integers 
if O is fermionic). Similarly, in the cases where O contains K types of matrices of 
different frequency Ui it can only connect states whose energy differs by ^^^niWi, 
where rii are integers whose absolute values are bounded by the number of matrices 
of each type appearing in O. 
As a result, ( p.8|) must have the form ( |3.6|) . Note that the argument based on ( |3.8|) applies 
not only to the thermal ensemble, but in fact to correlation functions in any density matrix 
(or pure state). 

To summarize, one finds that in free theory a real-time thermal two-point function 
always has a discrete spectrum and is quasi-periodic in the large limit. This implies 
that once one perturbs the theory away from thermal equilibrium, the system never falls 
back and keeps oscillating. This is not surprising since the system is free and there is no 
interaction to thermalize any disturbance. Note that this is distinctly different from the 
behavior (|2.15 ) found at strong coupling. In particular, this implies that the bulk descrip- 



tion of the high temperature phase in free theory looks nothing like a black hole. Also 
note that the story here is very different from that of the orbifold CFT in the AdSa / CFT2 
correspondence. There the mass gap in free theory goes to zero in the large N limit in the 
long string sector P^]. As a result, one finds that free theory correlation functions in the 



long string sector do resemble those from a BTZ black hole . 

3.2. Perturbation theory 

In this subsection we use a simple example ( p.3|) for illustration. The general features 
discussed below apply to generic theories in (|2.1J ) including A/" = 4 SYM. 
In perturbation theory Ge{t) can be expanded in terms of A as 



Ge{t) = Y,\-G^^\t) (3.9) 

T1 = 



where is the free theory result. We will be only interested in the connected part of 
Ge{t). Higher order corrections are obtained by expanding -f '^^^ in the path integral 



with V given by the quartic term in (|2.3|) . More explicitly, a typical contribution to G^ (r 
in ( 3.9|) has the form 



-11 



/ rfri---/ rfr„ (O(r)0(0)y(ri)---y(r^))^o (3.10) 
Jo Jo 
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The free theory correlation function inside the integrals in ( p.lO|) can be computed by first 
using Wick contraction ( |3.2| ) and then doing the U integral. The general structure of ( |3.10| ) 
can be summarized as follows: 

1. The planar diagram contribution to G^^\t) scales like A^*^, while diagrams of other 
topologies give higher order l/N"^ corrections. The number of planar diagrams 
grows like a power in n, i.e. is bounded by with C some finite constant ||34|| . 

2. The r-integrations are over a compact segment and are all well defined. A typical 
term in (|3.10 ) after the integration has the structure 

gif{(3ye''^<^- (3.11) 

where I and k are integers. I can take values from to n, while k from —2n — A to 
2n + A where A is the dimension of O in free theory. 

Analytically continuing ( |3.9D to Lorentzian time by taking r = zt + e, we find 

oo 

G+(t,A) = 5^A"Gf)(t) (3.12) 

n=0 

where typical terms in G^_^^ (t) have the t-dependence of the form 

g\2\fi)t^e'^"'' (3.13) 
with the range of / and k given after equation ( |3.11| ). After Fourier transforming to 
frequency space we find that at each order in the perturbative expansion G^^\u) (and 
thus the spectral density function p{uj)) consists of sums of terms of the form 

47)5«(^-fc^o) (3.14) 
where the superscript / denotes the number of derivatives. 

If the theory has more than one fundamental frequencies, since the interaction vertices 
are traces of a finite number of matrices they only connect states with definite energy 
differences. More and more frequencies will appear in the spectrum of a correlation function 
as we go to higher and higher orders in the perturbative expansion. The increase in the 
number of frequencies is exponential in the order of the expansion but at any fixed finite 
order no matter how large the spectrum of the correlation functions is discrete. 

One origin of terms in ( p.l3|) is the shifting of frequency from the free theory value. 
For example, suppose the free theory frequency is shifted to u = ujq + \ui + ■ ■ ■, one would 
get terms of the form ( p.l3| ) when expanding the exponential e*^* in A. One can in principle 
improve the perturbation theory by resumming such contributions using Dyson's equations. 
However, there appears no systematic way of doing this for a composite operator (|3.4|). 
In Appendix B.2, we prove that real-time correlation functions of fundamental modes Ma 
again have a discrete spectrum in the improved perturbative expansion. 
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4. Break down of Planar perturbation theory 

It is well known that at zero temperature the planar expansion of a matrix quantum 



mechanics has a finite radius of convergence in the A-plane (see e.g. for a recent 
discussion and earlier references). If this persists at finite temperature, properties of the 
theory at zero coupling or in perturbation theory should hold at least for the coupling 
constant being sufficiently small. In particular, from our discussion of last section, one 
would conclude that real-time correlation functions for generic gauge invariant operators 
should be quasi-periodic and an arrow of time does not emerge at small 't Hooft coupling. 
In this section, we argue that the planar perturbative expansion in fact breaks down for 
real-time correlation functions and thus perturbation theory cannot be used to understand 
the long-time behavior of the system at any nonzero coupling. 

From our discussion in section 3.2, we expect the Euclidean correlation function ( |3.9| ) 
should have a finite radius of convergence for any given r G (0, (3). After analytic continu- 
ation to real time, the convergence of the expansion in Euclidean time implies that ( |3.12| ) 
should have a finite radius Xc{t) of convergence for any given t. However, it does not tell 
how Xc{t) changes with t in the limit t —>■ oo. In this section we argue that the radius of 
convergence goes to zero in the large t limit. Note that the convergence of the perturbative 
expansion depends crucially on how g^^^^ in ( p.l4|) fall off with n. We will argue below that 
the falloff is slow enough that perturbation theory breaks down in the long time limit. In 
frequency space, one finds that n-th order term in the expansion grows like nli^l. 

We will again use ( |2.3| ) as an illustration. The argument generalizes immediately to 
generic systems in ( |2.1| ). For simplicity, we will consider the high temperature limit ( |3.4| ) 
in which we can replace U in ( |3.2| ) by the identity matrix, e.g. 

Miij(r) Mifcz(O) = — firj^^V - 171(3; ujQ)diidkj = ^^u^kj 9E{r;uj) (4.1) 



where 



with 



gEir; ^) = ^ (e-""(l + f{oo)) + e'^V(c^)) , r e (0, (3) (4.2) 



n^) = • (4-3) 



Note that in frequency space the relation between real-time and Euchdean correlation func- 
tions is not simple, since Euclidean correlation functions are only defined at discrete imaginary 
frequencies. 
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Note that outside the range in ([4.2|), QEij) is periodic. 

For our purpose it is enough to examine the Wightman function for Mi, 



D+{t) = (e-^^Mi(t)Mi(O)) . (4.4) 

An exactly paraUel argument to that of the last section leads to the expansion 

oo 

I)+(t,A) = 5^A-I)f)(t) (4.5) 

ri=0 

where typical terms in D^"* {t) have the t-dependence of the form 

4"^(/9)t^e^^'"°* (4.6) 

The convergence of series depends on how d!)^^ fall off with n. For our purpose it is enough 
to concentrate on the term with the highest power t in each order, i.e. the coefficients of 
t"^ with given k. More explicitly, we will look at a term of the form 



D+{t, A) = Df{t) J2 CnX'^t^ + ■■■ (4.7) 

n=0 



where D^"^ is the free theory expression. 



As before we will first compute ( f4.7|) in Euclidean time and then perform an analytic 
continuation. Calculating 

Cn explicitly at each loop order for all n is of course impractical. Our strategy is as 
follows. We will identify a family (in fact infinite families as we will see below) of planar 
Feynman diagrams of increasing loop order and show that their contribution to falls 
off like a power in n. Barring any unforseen magical cancellation0, this would imply that 
the perturbation series (|3.12[) has a zero radius of convergence in the t ^ oo limit. The 



simplest set of diagrams which meet our purpose are given by: 



Note that since we are in the high temperature phase, in which supersymmetry is badly 
broken, there is no obvious reason for suspecting such magical cancelations. 
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Fig. 1: A family of diagrams which indicates that the perturbation theory break 
down in the long time limit. Black and red lines denote propagators of Mi and M2 
respectively. 



These graphs appear at orders di = 2,d2 = 8, ds — 26, ■ ■ ■ of perturbation theory 
where 

d, = 3d,-i + 2 = 3^ - 1, z = 1, 2 ■ ■ ■ . (4.8) 

We denote the contribution of each diagram by ri(r). For our purpose, it is not necessary 
to compute the fuU graph. We wiU only need to calculate the term in each graph with the 
highest power of r, i.e. the term proportional to t'^\ Also note that in each diagram, the 
symmetry factor is exactly 1. Let us start with Fi, which is given by 

Fi(r-r') = A^/ dTidT2 gE{r - Ti;uo)gE{ri - r2;uo)gEiT2 - r';uJo) (4.9) 
^0 

Note the identity 



9Ur;^o) = J^f M (^e^"°^7E(r;cc;o) + J^9Eir;3uo)j (4.10) 

Now plug (|4.10| ) into ( [1.9|) . It is easy to convince oneself that the term proportional to 
(7£;(r; Scuq) in ( fl.lOp will not generate a term proportional to and we will ignore it. The 
contribution of the term proportional to 5'e(^; '^0) can be found by noting the identity^ 

/ dTidT2 gE{r - n; u;o)fi's(ri - T2]UJo)gEiT2 - r'; wq) = 77 .0 2\2 3e{t - t';ujo) (4.11) 
Jo ^ yo'^o) 



If two matrices have different frequencies in the product gE{T\u)o)g%{T\u)\) there is also a 
term proportional to gE{T;LUo) with coefficient ■(2^~)t/('^i)(1 + /(c<ji)) and the rest of the analysis 
follows with minor changes. 
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/ _ f\2 

The right hand side of ( [4 .111) contains a piece | (2ljo)^ dEi^ — t') plus parts with smaUer 
powers of r — t' . Thus the term in (WM) proportional to (r — r')^ is given by 



ri(r - r') = ^(r - r'fgE{r - r') + ■ ■ ■ (4.12) 

where 

3/(1 + /) , . 
« = ^^^, / = /M (4.13) 

The term proportional to r*^' for higher order diagrams Ti{T) can now be obtained by 
iterating the above procedure. A useful identity is 

/ dTidT2 gE{r - ri]UJo)gE{ri - r2; cuq) (n - r2)'^ 9E{r2 - r'; Uq) 
"^0 (4.14) 

(2^0)"^ (n + 2)(n + l) 
where we kept only the term with the highest power of r — r', as lower power terms will 
not contribute to the terms in which we are interested. We find that the term proportional 
to r'^" in ri(r) is given by 

r,(r) = F,X''^r''^gE{r; ujo) + ■ ■ ■ (4.15) 
where satisfy the recursive relation 

F.+i=F^——^ TT- (4.16) 

di+i{di+i - 1) 

Thus Fi can be written as 

F, = a^A, (4.17) 

with 



fc=0 



Ai in the large i limit can be easily estimated and we find 

Ai ^ e~^'^\ i > 1 . 



Summing all our diagrams together and analytically continuing to Lorentzian time 
with r = it + e, we find thatH 

j:m)^D^:\t)f^{-iy(^y + ■■■ (4.19) 



Since we are only interested in the asymptotic behavior of the sum for large i we have replaced 
Fi by its asymptotic value. 
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with he given by 



hr 



3 

62 



62 



(4.20) 



(4.21) 



v/3/(l + /) 

Equation (|4.19|) implies that the radius of convergence in A is given by 

Ac(t) - I 

which goes to zero as t ^ oo. 

It is also instructive to repeat the computation of fig. |l| in frequency space using 
the real-time method. The calculation is straightforward and we will only summarize the 
result. One finds that the contribution of Ti to the Feynman function Df{oj) grows like di\. 
Thus one expects that the perturbative expansion in frequency space is not well defined 
for any frequency. Note that the non-analyticity in frequency space can be expected since 
in going to frequency space one has to integrate the full real time-axis and the Fourier 
transform is sensitive to the long time behavior. Also note that the n! factorial behavior 
in perturbation theory often implies an essential singularity at A = (see also below). 

We conclude this section with some remarks: 

1. In the zero-temperature limit he —>■ oo and the set of terms in ( [4.19|) all go to zero. 

2. To simplify our discussion, we have only considered diagrams in fig. There are 
in fact many other diagrams of similar type contributing at other orders in A. For 
example, by including those in fig. 2, one can get contributions for all even orders in 
A rather than only p.8| ). The qualitative conclusion we reached above is not affected 
by including themii. 





Fig. 2: By including the diagrams on the left with all possible i,j,k > we can 
get contribution at every even order of A instead of ( [4.81 ). To denotes a single 
propagator. Diagrams on the right can also contribute to the odd orders if ( p.3[) 
contains additional interactions of the form txA^B^ . 



There are also potentially an infinite number of other sets of diagrams which can lead to the 
behavior ( ^^ ), e.g. one can replace Fi by any diagram whose highest power in t is the same as 
the order of perturbation and then iterates. 
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3. By taking in consideration the diagrams on the left of fig. ^ the sum in ( [4.19|) is ex- 
tended ii to aU even powers of Xt and is osciUating therefore the singularities in At 
should lie on the imaginary axis. Let us suppose that for a given A, -D+(t) has a sin- 
gularity mt at Qc/X with Qc lying in the upper half plane@. Now Fourier transforming 
Dj^it) we find that 

D+{uj) = / dte'^^D+{t) (4.22) 



The presence of Qc/X and q*/X implies D^{uj) contains a term of the form for ci; > 

D+(w)~e^'^^ (4.23) 



Thus D^(u) contains an essential singularity at A = 0. 
4. The n! behavior in perturbative expansion in frequency space (say in the computation 
of Df{uj)) arises from a single class of Feynman diagrams. This is reminiscent of 
renormalons in field theories []35| . |3^ , |37| . In particular, when Borel resumming the 



divergent series, depending on whether u is greater or smaller than wq? the singularities 
on the Borel plane can appear on the positive or negative real axisil, also reminiscent 
of the IR and UV renormalons. 

3 

5. Note that in the limit T — * oo, he in ( [4.20|) scales with T as /ic ~ -5^, i-e. 

Ac(t) ~ ^ (4.24) 
For fixed A, we expect a singularity for -D+(t) at 

t - ^ (4.25) 

Note that the right hand side of ( ^.25| ) is reminiscent of the magnetic mass scale for a 
Yang-Mills theory (see e.g. [0). However, in our matrix quantum mechanics, there 
is no infrared divergence and it is not clear whether there is a connection. 

6. The discussion can be straightforwardly applied to a generic theory in (p.lD with 
cubic and quartic couplings. In fact the argument also applies to a single anharmonic 



the value of he also changes 

Note that Qc/X must also be a singularity of D-^{t). 



Since we only have contributions to even order in A, we cannot make a conclusion from our 
discussion so far. 



19 



oscillator at finite temperature, even though in that case one does not expect the 
perturbative expansion to converge anywayii. Similarly, the argument also applies 
to a single-matrix quantum mechanics if one does not impose the singlet condition. 
When imposing the singlet condition, the matrix U in equation cannot be set 
to 1 and our argument does not apply. Similarly our argument does not apply to 
in the low energy sector, in which U is always important. Indeed using the 



results of [pl| , |39[| , one can show that to leading order in the large limit, correlation 
functions at finite temperature can be written in terms of those at zero temperature 
and we do expect that the perturbation theory has a finite radius of convergence. 
7. Our argument indicates that perturbation theory breaks down in the long time limit 
for a generic theory in (|2.1j ). However, for any specific theory (say A/" = 4 SYM theory) 
we cannot rule out magical cancelations which could in principle make the coefficients 
of n-th order term much smaller than indicated by the diagrams we find. If magical 
cancelations do occur in some theory, that would also be extremely interesting since 
it indicates some hitherto unknown hidden structure!!. 



5. Physical explanation for the breakdown of planar expansion 

In this section we give an alternative argument for the breakdown of perturbation 
theory in the long time limit, which complements that of last section. The discussion 
below should apply to a generic theory in (|2.1| ) . For definiteness we use J\f = 4 SYM as an 
illustration example. 

We first set up some notations. We write the full Hamiltonian as 

H = Ho + V{X) (5.1) 

with Hq the Hamiltonian of the free theory and V the interaction. We denote a free 
theory energy eigenstate by \a) with energy £„. |0) is the (unique) free theory vacuum. 

See Appendix D for further elaborations on the example of a single anharmonic oscillator and 
a discussion on the differences between the single anharmonic oscillator and the matrix systems 
under consideration. 

Since we are working at a finite temperature, supersymmetry alone should not be sufficient 
for the cancelations. 
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The energy eigenstates of the interacting theory H are denoted by |z) with energy Ei. |0) 
is the interacting theory vacuum. We can expand 



with Cia satisfying 



) = J2c,a\a) (5.2) 



a 



J]|Cza|' = J]|C.„P = 1. (5.3) 



8~0(n2) 



e~0(i) 

Fig. 3: The energy spectrum of free A/" = 4 SYM on is quantized. Typical 
degeneracy for an energy level e ~ 0{1) is of order 0(1). Typical degeneracy for a 
level of energy e ~ 0{N^) is of order e°(^ ). 



We first recaU some relevant features of the free theory energy spectrum of A/" = 4 
SYM on S^. Since in ( |2.1j ) are all integer or half- integer multiples ol uq = the free 
theory energy spectrum is quantized in units ^uq- Typical energy levels are degenerate. 
The degeneracy is of 0(1) in the low energy sector and of order e*^*-^ in the high energy 
sector. The exponentially large degeneracy in the high energy sector can be seen as follows. 
From ( |2.7| ) the density of states Oo(e) in the high energy sector is of order e^^^ \ Since the 
energy levels are equally spaced with spacings order 0(1), it must be that typical energy 
levels have a degeneracy of order e^^^ \ Alternatively, the number of ways to construct 
a state of energy of order 0{N'^) from 0{N'^) oscillators of frequency of 0(1) is clearly 
exponentially large in A^^. 

When a theory contains n > 1 incommensurate fundamental frequencies, the free 
theory spectrum at energies of order 0{N'^) will have level spacings of order 0{N~^^'^~^'>). 
Since the density of states are of order e"^*^^ \ the degeneracy of a typical state is again 
of order e*^*^^ as in the case of one fundamental frequency. That the level spacings go 
to zero as a negative power in A^ does not change our conclusion of previous sections 
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regarding the thermalization in free theory or at any finite order in perturbation theory, 
since as emphasized there a smaU operator can only connect states whose energy differences 
are of order 0{N^). Therefore such smaU level spacings cannot be accessed dynamically. 
The restriction above to a finite number of incommensurate fundamental frequencies is 
not essential. The conclusion applies to any theories in which the number of fundamental 
frequencies increases only as a power of the frequency. For these theories, the thermal 
ensemble is dominated by states built from oscillators whose frequencies are smaller than 
or of order of Thus the effective number of fundamental frequencies is finite. 

Now let us turn on the interaction V{X) ( p.2| ) with a tiny but nonzero A. We will focus 
on the high energy sector. Given that free theory energy levels are highly degenerate, one 
would like to apply degenerate perturbation, say to diagonalize F in a degenerate subspace 
of energy E ^ nN'^ and of dimension e^^^ \ For this purpose we need to choose a basis 
in the degenerate subspace. This is a rather complicated question, due to difficulties in 
imposing singlet conditionJi^. However, when fi is sufficiently large we expect the singlet 
condition not to play an important roleiil. So to simplify our discussion we will ignore the 
singlet condition below. A convenient orthonormal basis of energy eigenstates for Hq are 
then monomials of various oscillators (appropriately normalized), i.e. 

N 

nn • ^^-4) 

a i,j = l 

In the basis ( ^.4|) , if the full theory is not integrable, V can be effectively treated as 
an (extremely) sparse random matrixll. The sparseness is due to that each term in V can 
connect a given monomial state to at most A^'^ other states, where k is an 0(1) numbeJii. 
Randomness has to do with the large dimension of the subspace and to the fact that there 
is no preferred ordering for the states within the same subspace. Diagonalizing V, we thus 
expect, from general features of a sparse random matrix (see Appendix C for a summary). 



The trace relations are important for states of such energies. 

As remarked earlier, in the high temperature limit the saddle point for U (in (|3.2[ )) approaches 
the identity matrix. 

"^^ Here we restrict V to & single energy level. When including all energy levels V is banded 
and sparse. The banded structure is due to energy conservation. 

This is because each term in y is a monomial of a few matrices. 
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1. The degeneracy of the free theory wiU genericaUy be brokenll. 

2. A number of states of order e^'^^ •* will mix under the perturbation. 

3. The typical level spacing between energy levels should be proportional to the inverse of 
the density of states and is thus exponentially small, of order e~'^^^ \ Note that this 
is exponentially smaller than the level spacings in a free theory with a finite number 
of incommensurate fundamental frequencies. 

The story is in fact a little more intricate. We expect the degenerate perturbation to 
be a good guide if the spread of energy eigenvalues after diagonalizing V in a subspace is 
smaller than the spacings between nearby energy levels. The spread F of eigenvalues of V 
can be estimated by (see Appendix C) 

r'-rl = Y,\{a\V\b)\'^0{N') (5.5) 

for any nonzero A, where the sum restricts to a degenerate subspace. Note that (|5.5| ) only 
depends on that F is a single trace operator and does not depend on the specific structure 
of it. That r ~ 0{N) implies that it is not sufficient to diagonalize V within a degenerate 
subspace. It appears more appropriate to diagonalizeii it in a subspace with energy spread 
of order 0{N). Thus in addition, we expect that: 

4. an interacting theory eigenstate \i) is strongly coupled to free theory states \a) within 
an energy shell of order 0{N). 

This statement will be justified in the next section from a somewhat different perspective. 
That r ~ 0{N) for any nonzero A in the 't Hooft limit indicates a tiny A may not really 
be considered small perturbation after all. 

Various features discussed above when turning on a small A are clearly non- 
perturbative in nature. However, it may be hard to probe them directly using Euclidean 
space observables like partition functions and Euclidean correlation functions. These ob- 
servables probe only average behaviors within an energy difference range of order 0(T) or 

For Yang-Mills theories on , there are remaining degeneracies associated with the isometry 
group SO {4) of S^. Except when one considers the sectors with very large angular momenta on 
S^, typical representations of 50(4) are rather small and should not affect our general argument. 

This statement is of course only heuristic since there is no sharp criterion to decide what 
should be the precise size of the subspace. However, we expect the A'^ scaling should be robust. 
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larger and thus may not be sensitive to the changes in level spacings at smaller scalesll§. In 



contrast, real-time correlation functions are much more sensitive. For example, consider 
the Lehmann spectral decomposition of G+{t), i.e. 



where we have inserted complete sets of energy eigenstates \i) of the interacting theory. 
From ( [5.6| ), it is clear that G+{t) can in principle probe any small energy differences, 
provided one takes t to be large enough. This explains the breakdown of perturbation 
theory in the long time limit observed in G+{t). At large A^, the A — > and t —>■ oo limits 
do not commute. 

In this and the last sections we have presented two lines of largely independent argu- 
ments that suggest that planar perturbation theory breaks down in the long time limit. 
The first argument (last section) is based on an honest Feynman diagram calculation, which 
establishes the breakdown of perturbation theory, but does not tell us directly whether the 
spectral functions are continuous or not. The second argument (this section), based on 
the energy spectrum, is more heuristic, but implies that the spectral functions are con- 
tinuous in the large limit. There are reasons to believe that the two arguments should 
be closely related. As stressed earlier, the breakdown of perturbation theory from the 
class of Feynman diagrams considered in the last section only happens at a sufficiently 
high temperatures, at which the thermal ensemble is dominated by states of energy of or- 
der 0{N'^) and the energy spectrum becomes quasi-continuou^. Nevertheless, a precise 
relation between the two arguments is not clear at this point, as they use very different 
languages (one Feynman diagrams and the other energy levels). It would be very desirable 
to find a direct connection and to have an understanding of the time scale ( [4.25|) from the 
point of view of the energy levels. 

Of course if one is able to compute Euclidean observables exactly, one should be able to 
extract all the interesting physics. After all, real-time observables can be obtained from Euclidean 
ones by analytic continuation. It is just often the case that real-time physics is encoded in a very 
subtle way in Euclidean observables. 

Applied to the theory at zero temperature or at a temperature below the deconfinement 
temperature, the class of diagrams gives a convergent contribution. At such a temperature, the 
thermal ensemble is dominated by states of energy of order 0(A^°), which are not quasi-continuous 
when turning on interactions. 



G+(t) = -IV(e-^^O(t)O(0)) 




(5.6) 
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6. A statistical approach 

The argument of section 4 shows that the planar perturbation theory breaks down in 
the large time limit, but it does not tell us what the long time behavior is. Non-perturbative 
tools are needed to understand the long time behavior of real time correlation functions in 
the large limit. Here we develop a statistical approach, taking advantage of the extremely 
large density of states in the high energy sector. In this section we outline the main idea and 
the results, leaving detailed calculations to various appendices. The statistical approach 
enables us to derive some qualitative features satisfied by the Wightman function for a 
generic operator at finite temperature, including that it has a continuous spectral density 
function and should decay to zero in the long time limit. The features we find here are also 
shared by the Wightman function at strong coupling found from supergravity analysis. 

Our starting point is the Lehmann spectral decomposition of G'+(t) ( p.6|) , 

G+w = ^E^"'^^^^^^'"''^-^v.. (6.1) 

where 

p., = |(z|O(0)|j)|2 = |O,,f (6.2) 

In momentum space 

G+{u) = I ^e-^^»5(a; + E^)p,^ . (6.3) 

Matrix elements Oij can in turn be expressed in terms of those of free theory using (c^a 
was introduced in (|5.2| )) 

O,, = (z|O(0)b-) = J2clc,, {a\0\b) = J^clcjbOat (6.4) 

a,b a,b 

where we have inserted complete sets of free theory states and Oat = (oj|C'(0)|6). 

Since for sufficiently high temperature, the sums in ( |6.3|) and ( |6.4|) are peaked at 
an energy with an extremely large density of states, one should be able to obtain the 
qualitative behavior of pij and G^{u!) from statistical properties of Oab and Cia- As 
discussed in the last section, in the interacting theory, we expect typical level spacings scale 
with N like e~'^^^ \ In the large N limit, Ei can be considered as taking continuous values. 
Note that this by itself does not imply that Gj^{uj) has a continuous spectral decomposition, 
since it is possible that pij only has support for states with finite energy differences. We 
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argue below that pij has nonzero support between states with any Ei — Ej G (—00,00), 
which is independent of A^, and thus Gj^{uj) does have a continuous spectrum. 

Let us first look at the statistical behavior of Cia- For this purpose, consider the 
following density functions 

Pa{E)=Y,\c^a?5{E-E,) (6.5) 
i 

X^(e) =$^|cza|'5(e-e„) (6.6) 

a 

Pa{E), first introduced by Wigner is also called the local spectral density function or 
strength function in the literaturell. Using normalization properties of c^a, one finds that 



jdEpa{E) = l, Jde 



X^(e) = 1 (6.7) 



Pa{E) can be considered as the distribution of interacting theory eigenstates of energy E 
coupling to a free theory state \a). Similarly, Xi(e) gives the distribution of free theory 
states of energy e coupling to an exact eigenstate The mean and the variances of the 
two distributions are given by 



Ea = / dEEpa{E) = {a\H\a) {Ql 



ra = Tl = JdE{E- Eaf Pa{E) = I 



(6.9) 



de€Xr{e) = E,-{i\V\i) (6.10) 



(6.11) 



= A2 = jdeie- e,f x{E,, e) = J] | {t\V\j) 

Ea and Fa give the center and the spread of interacting theory energy eigenstates coupling 
to a free state \a). Similarly, and give the center and the spread of free theory 
states coupling to an interacting theory energy eigenstate Fa can be considered as a 
measure of correlation among energy levels of the interacting theory (since states whose 
energies differ by Fa could couple to the same free theory state and are thus correlated). 
Ai characterizes the range of free theory states which are mixed by perturbation. Note 
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These density functions have been frequently used in quantum chaos literature, see e.g. |41] 
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that the heuristic discussion after equation (5^) imphes that ~ 0{N), which we wiU 
confirm below using a different method. 

Individual energy eigenstates are rather hard to work with. We will consider mi- 
crocanonical averages of (|6.5|) and After all, for ( |0| ) and (|6.3|) we only need the 

behavior of pij averaged over states of similar energies. We will denote the average H of 
Xi(e) over interacting theory states of energy E by X£;(e) and similarly the average of Pa{E) 
over free theory states |a) of similar energy e by Pe{E). Since the averages involve a huge 
number of states and the large limit is like a thermodynamic limit in the high energy 
sector, we will assume that Xe{,^) is a smooth slow functionE§ of i.e. it depends on E 
only through E/N'^. Similarly Pe{E) is assumed to depend on e only through e/N"^. The 
center and variance of X£'(e) and Pe{E) will be denoted by e{E), S(E') = A^(i?), E{e), 
and cr(e) = r^(e) respectively^. These quantities should also be slow functions of i? or e 
as they inherit the property from Xe{^) and Pe{E). In the Appendix E we estimate these 
quantities and find that 

-e{E)=N^g{\E/N'') 
^{E) = N'^h{X,E/N^) 

_ (6.13) 

E{e)=N''g{X,e/N^) 

We emphasize that the large A^ scalings above only depend that V is given by A^ times 
single trace operators. Given that the underlying theory is not integrable and the extremely 
large number of states, we will thus approximate Cia for fixed z as a random unit vector 
which centers at li with a spread of order ~ 0{N). 

Now we turn to the statistical properties of Oab- Our earlier discussion for V in the 
free state basis ( p.4| ) can be carried over to any operator O of dimension 0(1). Thus Oab 



More explicitly, the average can be written as 

Xi^(e) = ^ Yl Me) (6.12) 

Eie{E-S,E+S) 

where 5 is small enough that Q.{E) does not vary significantly in the range [E — 5,E + 5). 

Note that a function f{E) is considered a slow function if it can be written in a form 
f{E) = N°'g{E/N'^), where g{x) is a function independent of A'^. 

which can also be obtained by the average of various quantities (|6.8D -( |6n ) to leading order 
in large N. 
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can be considered as an sparse banded random matrix. The matrix is banded since from 
energy conservation O can only connect states whose energy difference is smaUer than 
the dimension of O. Note that even though Oab is sparse, for each row (or column), the 
number of nonzero entries grows with as a power. 

To summarize, we will assume the following statistical properties for Cia and Oab'- 

1. For a given i, Cia is a random unit vector with support inside an energy shell of width 
0{N). In particular, the Cia satisfy the same distribution for \a) of the same energy. 

2. Oab is banded sparse random matrix, with the number of nonzero entries growing 
with as a power. 

Now consider any two states \i) and \j), with energies Ei and Ej respectively, for which 
uj = Ei — Ej ~ 0{1). One finds that — ~ 0(1) and the energy shells of the two states 
overlap significantly. Given that the number of nonzero entries in a row or column of 
Oab grows with as a power and that each element of Cia satisfies the same distribution, 
one concludes from ( |6.4| ) that Oij should have support for any iv = Ei — Ej ~ 0{1) and 
G+{uj) has a continuous spectrum for iv G (— cxd, +oo). Note that the fact that A ~ 0{N) 
is crucial for having a continuous spectrum uj G (— oo,+c)o). Suppose A ~ 0(1), the 
spectrum cannot extend to ±00 due to energy conservation. 

One can further work out more detailed properties of pij. Leaving the detailed cal- 
culation in various appendices, we find that (after averaging pij over states of similar 
energies) 

PE,E, = ^^(^; E) = e-5(^)A(a;; E) (6.14) 

where Q{E) and S{E) = logO(i?) are the density of states and entropy of the interacting 
theory respectively and 

El -f- E2 

E = , u! = E-\ — Eo ■ 

2 



Equation ( p. 14] ) is derived in Appendix G along with properties of A{uj; E) stated below. 



Some useful formulas used in the derivation are collected in Appendix F. A{uj; E) can be 
expressed in terms of an integral of Xe{^) and e{E) (see equations (G.3) and (G.8)) and 
satisfies the following properties: 

1. A{u!; E) is an even function of i.e. 

A{-uj;E) = A{uj;E) (6.15) 

2. As a; — > 00 

A(u;;i?)oce-^^(^)H, m = ^ (6-16) 
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3. A{uj, E) is integrable along the real axis and can at most have integrable singularities 
of the form 

A{oo- E) oc ^— , a, < 1 . (6.17) 

\UJ — UJs\ " 

4. Ae{u)) depends on E only through E/N^, i.e. it can be written as 

A{uj;E) = A{u:;ij), ^ = ^ (6-18) 

and A is a function independent of A^. 
Note that property 2 implies that in the large limit, PE1E2 ~ Ei — E2 ^ N"" with 
a > 0. 

The expression for ^+(0;) in momentum space can now be obtained by plugging ( |6.14| ) 
into ( |673| ) and using a saddle point approximation. We find that 

G+{u;) = ^ I di?e-^^e^(^)+^(^+-)e-^(^+-/2)A(^,i?/iV2) 
= e 2 A{u,np) 



(6.19) 



where 

_ E^ dS{E) 



(3 . (6.20) 



Note that since in the large limit, E can be treated as continuous and pij has support 
for any energy difference, it is appropriate to approximate the sum in (|6.3| ) by an integral. 
Also from the second line to the third line we have used that the quantity inside the 
bracket depends on E slowly and performed a saddle point approximation. We conclude 
this section with some remarks: 
1. Gj^ioj) has a continuous spectrum with ui G (—00, +00) in the large limit (note that 
u does not scale with A^). 



2. Since A{uj^ p) can at most have integrable singularities of the form (|6.17|) on the real 
axis, after a Fourier transform to coordinate space, Gj^ [t) must go to zero in the limit 
t — > 00. If A{lo; n) is a smooth function on the real axis, then G+{t) must decay 
exponentially with time. 

3. Considering the last line of (|6.19| ) definition for A{uj; p), for Af = 4 SYM on 

at strong coupling, the corresponding A{u!; fx) can be found by solving the Laplace 
equation for a scalar field in an AdS black hole geometry and be expressed in terms 
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of boundary values of renormalizable wave functions for the scalar field . In par- 
ticular, A{u; fx) found at strong coupling satisfy all the properties ( |6.15| )-( |6.18|) (it is 
a smooth function on the real axis). 
4. It should be possible to obtain an explicit expression for A{u!; fx) (and thus G+{u!)) 
using the expressions found in the appendices (e.g. equation (G.3)) if one can find 
the density functions ( |6.5| ) and ( |6.(i| ) for a sparse banded random matrix with varying 
density of states. While those for constant density of states have been discussed in the 
literature (see e.g. not much appears to be known for the non-constant density 

of states [pj. 



7. Discussions 

In this paper we first showed that in perturbation theory, real-time correlation func- 
tions in the high temperature phase of (2.1) have a discrete spectrum and the system does 



not thermalize when perturbed away from thermal equilibrium. We then argued that the 
perturbative expansions for real-time correlation functions break down in the long time 
limit. The breakdown of perturbation theory indicates that at large the A — and 
t ^ oo limits do not commute. The reason for the breakdown is that a wide energy range 
(of order 0{N)) of degenerate free theory energy eigenstates mix under the interaction. 
The level spacings in the energy spectrum of 0(1) in the free theory become e~^^^ \ As 
a result, real-time correlation functions develop a continuous spectrum for any nonzero A. 
The continuous spectrum was argued from a statistical approach developed in section 6, 
where we also show that real-time correlation functions should decay to zero as t ^ oo 
and the system becomes time irreversible. 

We should emphasize that our arguments in this paper are qualitative in nature and 
far from foolproof. For example, instead of being a random vector, Cia could have some 
structure (e.g. being very sparse) within the range of its spread, in which case our statistical 
argument will not be valid. 

It is also important to emphasize our results only apply to the high energy sector and 
in the low energy sector (or in the low temperature), there is no indication of breakdown of 
the planar expansion. In particular the results we describe here are not inconsistent with 
that the sector near the vacuum might be integrable in the large limit [^^. In fact the 
results of are consistent with that the interpolation between the free theory and 

strong coupling may be smooth in the low temperature phase. 
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Our results indicate that there is a large "phase transition" at A = 0, i.e. physical 
observables undergo qualitative changes in the limit A 0. The "phase transition" we 
find here is somewhat unusual, since it is not manifest in the Euclidean quantities like 
the partition function. The partition function appears to be smooth in the A ^ limit. 
The "phase transition" is in real-time correlation functions and their Fourier transforms. 
Real-time correlation functions decay to zero at large time at any finite A, while oscillatory 
for A = 0. In frequency space there is an essential singularity at A = 0. 

It would be interesting to understand whether one can continue the physics at small A 
to large A. If there is no further large "phase transition" in A, we expect that the analytic 
structure of various correlation functions observed at strong coupling should also be present 



at small A. Such structure include the signatures of black hole singularities [ p9| , |28[] and the 
bulk-cone singularities p4| |. 

Given that an arrow of time emerges for small A in the large limit, it is natural 
to ask what should be the string theory interpretation of the high temperature phase for 
A/" = 4 SYM on at weak coupling, or from the micro canonical point of view, what is 
the bulk interpretation for a generic state in the high energy sector. 

From the parameter relations in AdS/CFT, 

1 Gn 1 



s 



N i- ^ ,8 „2 



Gn = 11-9X. ■ (7.1) 



RS Ar2' 

one might conclude that at weak coupling A ^ 1, ^ i? , i.e. the string length Is is 
much bigger than the AdS curvature radius R. However, it seems unlikely one can give an 
invariant meaning to the statement. For example, even starting with a metric with i? <^ /g, 
one could perform a field redefinition of the form g^j^^, g^j^^, + a' R^^, + • ■ ■. In terms of new 
metric one then has R Ig. Thus it seems to us that even for A ^ 1, the corresponding 
bulk string theory should describe a spacetime of stringy scale, rather than sub-stringy 
scale. This is also expected from the gauge theory point of view. At weak coupling the 
only mass scale is the inverse radius of the sphere and there are no other lighter degrees of 
freedom. Thus the string scale has to be of the same order as that of the AdS curvature 
scale. 

Can one interpret the bulk configuration corresponding to the high temperature phase 
at weak coupling as a stringy black hole? It seems to us the answer is likely to be yes. Let 
us list the properties that the corresponding bulk configuration should satisfy as expected 
from gauge theory, assuming there is no further large "phase transition" between weak 
and strong couplings: 
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1. The bulk configuration should have an entropy and free energy of order 0{l/g1). 

2. The object absorbs all fundamental probes (since boundary correlation functions decay 
with time). 

3. The bulk geometry should have a horizon (since the boundary theory has a continuous 
spectrum) . 

4. The bulk configuration is likely to have singularities (since the signatures of the black 
hole singularities in gauge theory at strong coupling cannot disappear as the coupling 
is changed if there is no phase transition). 

5. A generic matter distribution will collapse into such a configuration (since in the 
boundary theory, a generic initial state will approach the thermal equilibrium) . 

6. Results in [ P^ , |2B[ ] indicate that the Euclidean time circle in the dual geometry for the 
theory in the high temperature phase should become contractibleEl. 

From the properties above, it seems appropriate to call it a stringy black hole. 

Finally let us mention that it is possible that a stronger version of equation ( |6.14| ) 
holds, i.e. for two generic states \j) in the high energy sector, 

P^J = -^M^■.E)■R,J, (7.2) 

with 

Ei -\- Ej 

E = , u) = Ei — Ej 

2 



and IZij a random matrix. Equation ( [7.2| ) is considered to be the hallmark of quantum 



chaos p5| , ^ il. Thus it is possible that A/" = 4 SYM is chaotic in the high energy sectoJll. 
Such a chaotic behavior, if it exists, might be related to the BKL behavior near a spacelike 
singularity [0. 
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Note that this alone cannot imply that the bulk geometry is a black hole since even at 
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It has also been argued in |46] that if (7.2) holds, then thermalization always occurs. 
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Appendix A. Parameter relations in AdS/CFT 

Af = 4 super Yang-Mills theory is a conformally invariant theory with two parameters: 
the rank of gauge group and the 't Hooft coupling A = gy^N. From the operator-state 
correspondence, physical states of the theory on can be obtained by acting with gauge 
invariant operators on the vacuum and their energies are given by the conformal dimensions 
of the corresponding operators. 

It was conjectured in that J\f = 4 SYM gives a nonperturbative description of type 
IIB superstring theory in AdS^ x 5*5. The AdS string theory also has two parameters: the 
ratio between string length Ig and the curvature radius R of AdS, and the ratio between 
the (lOd) planck length Ip and R. These ratios respectively characterize classical stringy 
corrections and quantum gravitational corrections beyond the classical supergravity. For 
small Ig/R and Ip/R, parameters of SYM theory and bulk string theory are related by^ 

- — , G'jv = L~£?,a . (A.l) 



yx' R^ N^'' 



The above relations indicate that the classical supergravity limit is given by the large N 
and large A limit of the SYM theory. In particular, a departure from the large N limit of 
the Yang-Mills theory corresponds to turning on quantum gravitational corrections in the 
AdS spacetime, while a departure from the large A limit (with N = 00) corresponds to 
turning on classical stringy corrections. 

AdS/CFT implies an isomorphism between the Hilbert space of the two theories. 
In particular, any bulk configuration with asymptotic AdS^ boundary conditions can be 
associated with a state (pure or density matrix) of the Yang-Mills theory. The mass M of 
the bulk configuration is related to the energy E in YM theory as 

E ~ MR . (A.2) 



"^^ We omit order one numerical constants. 
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Depending on how E scales with in the large limit, states of Yang-Mills theory are 
related to different objects in string theory in AdS. For example those whose E do not scale 
with (i.e. of order should correspond to fundamental string states. An object in 

AdS with a classical mass M satisfies 

GnM = fixed, Gn/R^ (A.3) 

From ([A.l| ) and ( [A .21 ) the corresponding state in YM theory should have E ~ 0{N'^). 



Appendix B. Self-energy in the real time formalism 

In this appendix we first review some basic properties of real-time correlation func- 
tions. We then prove that the spectral density functions of fundamental fields in (|2.1|) have 
a discrete spectrum after the resummation of the self-energy diagrams a la Dyson. 

B.l. Analytic properties of various real-time functions 

Various real-time thermal Wightman function for an operator O are defined by 

G+it) = ^Tr{e-f'^O{t)O{0))-C 
G_{t) = ^Tr {e-^^0{Q)0{t)) - C 

Grit) = e{t)G+{t) + e{-t)G-{t), (B.l) 
Gn{t) = ze{t)^Tv {e-^^[0{t),Om), 

GA{t) = -te{-t)^Tr{e-^^[0{t),Om) 

where Z is the partition function and C is a constant to be specified below. It is also 
convenient to introduce 

Gi2{t)=G+{t-i(3/2) (B.2) 

which can be obtained from G+(t) by an analytic continuation. 

By inserting complete sets of states in ( p.l| ), G'+(t) can be written as 

G^it) = ^Y.'-^^''^''''''^''' P^J (B.3) 
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where i,j sum over the physical states of the theory and pij = \ {i\O{0)\j) |^. Comparing 
( |B.3| ) and ( [B.l|) , C is chosen to be 

C = ^Y.'-"^'' P^^ (B.4) 

i 

Note that C is chosen so that the Fourier transform of Gj^{t) does not have a "contact" 
term proportional to 6{uj). Assuming the convergence of the sums is controlled by the 
exponentials, it follows from ( [B.3|) that is analytic in t within the range — /3 < 

Imt < 0. Similarly G-{t) is analytic for < Imt < /3 and 6*12 (t) for -| < Imt < |. 
Introducing the spectral density function 

then the Fourier transforms of flB.ll) can be written as 



2smh^ 



Gr{u) 



-00 
00 



2 

P(^0 (B.6) 
27r u — u' + it 

duj' p{i^') 



3 27r — cu' — ze 
Gf{uj) = GR{uj)+iG-{uj) 



Prom (IB. 61) we also have 



pM = -z(G'flH-G'AH) (B.7) 

We also note that the Euclidean correlation function in momentum space can be obtained 
from 

'Gr{zuji) l>0 27rZ 
Ge{uji) = { , ui = —, leZ (B.8) 



GAitcoi) / < ^ ' 



Some further remarks: 
1. From ([B3|)-([B:6D, 



p{-uj) = -p{u), G^2{-oo) = G12M, Gr{-uj) = Ga{oo) . (B.9) 

2. For a theory with a discrete spectrum, from ( |B.5| ), the spectral function p{u}) and 
Gj^{uj) are given by a sum of discrete delta functions supported on the real axis, while 
Gr{uj) is given by a discrete sum of poles along the real axis. 
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B.2. Self-energy in real-time formalism 

In this section we consider real-time correlation functions of fundamental fields Mq, 
in perturbation theory using the real-time formalism. We denote various quantities in 
( [B.l| ) with O = Met by D^-*, D^-* etc and will suppress superscript a from now on. We 
prove that the corresponding spectral density functions have a discrete spectrum after the 
resummation of the self-energy diagrams a la Dyson. For simplicity, we will consider the 
high temperature limit so that we can ignore the singlet projection (see (|3.4| )). 

In the real time formalism 13^] the degrees of freedom of the theory get doubled (see 
also [0). For each original field (type 1) in ( pTT|) one introduces an equivalent field (type 
2)E1 whose interaction vertices differ by a sign from the ones for fields of type 1. Vertices 
therefore do not mix the two different kind of fields but propagators do and are written as 
a 2 X 2 matrix. For example, in frequency space the propagator for Mq, (in the interacting 
theory) can be written as Dab{oj), a,b = 1,2 with each component given by 



13 

Dab can be diagonalized as 



(B.IO) 



with 



Am 



and 



jj ^ I cosh7 sinh7 \ ^^^^ ^ ' ^.^^^ ^ - /g 

\^sinh7 cosh7y'' ^ Ve^l^l - l' ^ Ve/^l^l - 1 ^ ' ' 

D,{u) = z / ^^Sp = J (B.13) 

The last expression in ( p.l3| ) implies that when analytically continued from the positive 
real axis, Dg{uj) cannot have singularities in the upper half w-plane. Similarly when 
analytically continued from the negative real axis, Dg{uj) cannot have singularities in the 
lower half w-plane. Note that Dg{ui) can have a discontinuity at Im{uj) = 0. If Dg{uj) does 
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In a path integral derivation these correspond to the fields whose time argument \s t — ia 



and we will take a — 2 • 
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turn out to be analytic on the real axis, then it can have singularities only on the real axis 
in the limit e ^ 0, which in turn implies that Dfi, Da and Dp can have singularities only 
on the real axis in the limit e ^ 0. 

We will now show that Dg{uj) obtained using the Dyson equation from any finite 
order computation of the self-energy is a rational function with singularities only on the 
real axis. This implies that the spectral function p consists of a sum of finite number of 
delta functions supported on the real axis. 

Note that the Dyson equation can be written as 

^ ^ -ili{u) (B.14) 



where 

ijj'^ — + It 

is the free theory expression and T^{uj) can be computed from perturbation theory as 
follows: (i) Compute 2x2 matrix Ilab{oj) from the sum of amputated IPI diagrams for 
the propagator in real time formalism; (ii) Diagonalizing Ilab{^) using ( p.l2|) , i.e. 

n„, = t/f^M ^\\u. (B.16) 
V n*(w) J ^ ^ 



That liah can be diagonalized using C/ is a consequence 
Now expanding Dg and 11 in power series of A 

Dg = Df) + AD(i) + A^D(^) + ■ • • 
^ Z 1 B.17 

n = An(i) + A2n(2) + . . . 

from equation ( p.l4| ) we have 

= ^ ^(o)(^n(2))^(o) ^ i^(o)(in(i))2;(o)(^n(i))j)(o)^ . . . 

(B.18) 

^From our discussion in section 3 (applied to fundamental fields), at any finite order in 
perturbation theory p{uj) consists of sums of terms of the form (|3.7|). Plugging such a p{uj) 
into ( p.l3| ) one finds that Dg^\uj) is a rational function and is analytic on the real axis 
at each order in the perturbative expansion (i.e. there is no discontinuity at Im{uj) = 0). 
Using ( |B.18| ) we find that 11*^ [u) must also be a rational function and analytic on real axis. 
This in turn implies that the resummed Dg{uj) found from (p.l4| ) is a rational function 
and analytic on real axis. We conclude that the singularities of Dg must lie on the real 
axis and there are only a finite number of them at any finite order in the computation of 
the self-energy H. From ( [B.13D the spectral density function must be a finite sum of delta 
functions supported on the real axis. 
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Appendix C. Energy spectrum and eigenvectors of sparse random matrices 

In this appendix we summarize features of eigenvalues and eigenvectors of a random 
sparse matrix found in p2| , |53| , pl| . Consider an M x M real symmetric matrix A whose 
elements Aij for i > j are independent identically distributed random variables with even 
probability distribution f{Aij). Let f{x) be of the following form: 

fix) = {l-a)6{x) + ah{x) (C.l) 

where < a < 1 and h{x) is even and not delta-function like at x = 0. Let the variance of 
h{x) be v'^. The parameter a measures the sparsity of the matrix: for each row or column 
of the matrix there will be on average aM = K elements which are different from zero. K 
is called the connectivity of the matrix. When iiT < 1 it is possible for eigenvectors to be 
localized in a subspace with dimension smaller than M. For i^T > 1 in the large M limit no 
such localization occurs and the matrix has to be diagonalized in the full M dimensional 
space. 

When K ^ 1, the density of states reduces to Wigner's semicircular law in an expan- 
sion in K~^: 

P{E) = ^v/4r2-i?2 (1 + 0{K-')) (C.2) 
Where E is the eigenvalue value and F is given by: 

= Kv^ {1 + 0{K-^)) (C.3) 

Notice that Kv^ is the average value of 

r? = $^|A.,f (C.4) 

over the rows or columns of the sparse matrix. The first correction to p{E) gives a change 
in the edge location, however there also are nonperturbative tails to the distribution which 
for E ^ r assume the form: 

(C.B) 

Their effect is to make the spectrum unbounded. 

Denote by T the orthogonal change of basis matrix which brings A to diagonal form 
for K ^ 1. T has a random uniform distribution over the group of orthogonal M x M 
matrices. Therefore the eigenvectors of A are a random orthonormal basis of the total 
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space which means that apart from correlations!! which are negligible in the large M limit 
their elements are independently distributed gaussian random variables with mean and 
variance j^. In particular the eigenvectors are completely delocalized. Therefore for large 
K the situation is similar to that for the Gaussian Orthogonal Ensemble (GOE). 



Appendix D. Single anharmonic oscillator 

It is clear that the argument presented in section 4 applies to the real time corre- 
lation functions of a single anharmonic oscillator at finite temperature (with changes of 
combinatorial factors) 

S^Jdt Qx^ - - Iax^) . (D.l) 
For example one can conclude that the perturbation theory for 

D+{t) = (x(t)x(O))^ (D.2) 

should diverge at a time scale ( [4.25| ) for T ^ t^o (we set ci^o = 1 in (p.lj )). Here we give 
an alternative derivation of this. Inserting complete sets of states in (|D.2| ) we find that 

D+{t) = Z-^Y,\ (^l^^lw) l^e-^^"-'*^^---^") (D.3) 

n,m 

where Z = ^-PE^ interacting theory eigenstates. If we are interested only 

in contributions of the form {\t)^ we get: 

D+{t) = Z-' Yl I (^l^l^)o |2e-/^^i°'-^*(^^^°'-^i°')-*^(^^''-^i'') + . . . (D.4) 

where quantities with index are computed in the free theory and XE^'^^ are the energy 
shifts at first order in perturbation theory. Equation ( p.4| ) can be evaluated as 

^ oo 

D+{t) = -Z-^ ^(n + i)[e-/5(-+i)-*(i+i(-+i)) + e-/3(n+|)+zt(i+A(„+i))j . . . 

n=0 

which can be summed to give 

(g/^ - l)e/^-**(i-#) (e/^ - l)e^*(i-i) , , 



which are due to normaUzation conditions. 

This section is motivated from a discussion with Steve Shenker. 
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In ( p.6|) there are double poles at 



t = ±z^ + fci^, kez. (D.7) 
A A 

If one resums the diagrams discussed in section 4, one would then get simple poles and the 
positions of the poles are further away from the real axis than those of (|D.6|) indicating 
that there are some positive contributions not captured by the class of Feynman diagrams. 

The reason for the behavior ( p.5| )-( p77|) can be attributed to the fact that the first 
order energy shift behaves as 

A(£;2i-i?(i))ocAn. (D.8) 

Thus when n is sufficiently large i.e. n ~ perturbation theory breaks down due to level 
crossing. Also note that the divergence of perturbation theory at t ~ ^ has nothing to do 
with the standard argument of the breakdown of perturbation theory by taking A —A. 
Indeed the behavior here is due to a single class of diagrams not to the n\ growth of the 
number of diagrams. 

We emphasize that while from the Feynman diagram point of view the discussion for 
anharmonic oscillators is almost identical to that for a matrix quantum mechanics (ex- 
cept that for matrix quantum mechanics one restricts to planar diagrams) , the underlying 
physics for the breakdown of perturbation theory appears to be very different: 
1. In the example of a single anharmonic oscillator, perturbation theory is asymptotic, 
i.e. the n-th order expansion contains n\ independent diagrams. In contrast, in the 
planar expansion of a matrix quantum mechanics, the number of Feynmann diagrams 
grows only like a power in n. The class of planar diagrams we identified gives rise to 
0{n\) contribution (in frequency space) at the n-th order. In the anharmonic oscillator 
example, given that the perturbative expansion is already divergent, one cannot really 
draw any clear conclusion from this class of diagrams. For instance, its contribution 
could be overwhelmed by those from n! other diagrams. In contrast, in the case of ma- 
trix quantum mechanics, the contribution from the particular class of diagrams makes 
an otherwise convergent perturbative expansion divergent. Given that the nature of 
perturbation theory is very different between the single anharmonic oscillator and the 
planar matrix quantum mechanics, one should be very careful in drawing any con- 
clusion when comparing them. In particular, the fact that the anharmonic oscillator 
has a discrete spectrum does not imply that in the matrix quantum mechanics case. 
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the divergence of the subclass of diagrams and hence the breakdown of perturbation 
theory are not related to a possible underlying quasi-continuous spectrum. 
2. As indicated earlier in this appendix for a single anharmonic oscillator, the divergent 
behavior of the class of Feynman diagrams considered in section 4 should have to do 
with with level mixing for states of energy 0(1/A). Applying the same technique to a 
matrix quantum mechanics, one again expects to relate the divergent behavior of the 
class of Feynman diagrams to the mixing of energy levels which dominate the thermal 
ensemble (i.e. with energy 0{N'^)). More explicitly, let us write ( p.4| ) for the matrix 
case as 

D+{t) = Zo-i 5^e-^<°' I {^\M\m), |2e-^*(^^.°'-^i"^)-*Mi^L^^-i^i^') + . . . (d.9) 

n m 

As we discussed in the main text, the sum over m in the above equation will involve 
an exponentially large number of states with free theory energies ranging over of order 
0{N). A naive estimate of —E^ '^ also gives order 0{N). Here unfortunately the 
story appears to be rather complicated and it appears it is not possible to extract a 
divergent time scale 1 / AT from (|D.9| ) . 

In summary, the Feynman diagram argument demonstrates the breakdown of pertur- 
bation theory, but does not tell us why or how it breaks down. It is certainly possible that 
completely different mixing behaviors in the energy levels may be reflected similarly by 
Feynman diagrams. In the anharmonic oscillator example discrete levels mix, while in the 
matrix quantum mechanics a quasi- continuous spectrum mixes. One must be very careful 
in extrapolating the results for an anharmonic oscillator to a matrix quantum mechanics. 

Appendix E. Estimate of various quantities 

We now estimate (|6lD-(|Jl|) after averaging them over states of similar energies. 
We will be interested in how these quantities scale with in the large N limit. An 
important property that we will assume below for these averaged quantities is that they 
are slow- varying functions of e or i?. In the large limit, we can then estimate them using 
the corresponding thermal averages, which can in turn be expressed in terms of various 
correlation functions at finite temperature. For example, the thermal average of is 

= |5^e-^^^E, = 1 /■ dEe-^^n{E)J:{E) (E.l) 

i 
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where is the microcanonical average and fl{E) = e^^^^ the density of states. Since 

E(i?) is a slow- varying function of E, we can perform a saddle point approximation of the 
last expression, yielding 

E{E) ^ S(/3e) (1 + 0(l/iV2)) (E.2) 

with determined by = Pe- Using the last equality of ( |6.11D we can write 

as 

= {v{o)vm^ 

where {V{0)V{0))p denotes the connected Wightman function as defined by ( p.3|) . From 
the standard large scaling argument ( [E.3|) is of order 0{N'^) (recall that we include a 
factor of in the definition ofV). Thus unless ( [E.3|) is zero at leading order we conclude 
that can be written in a form 

E{E) = N^h{X, E/N^) (E.4) 

where /i(A, fj.) is a function independent of A^. An exactly parallel argument can be applied 
to a(e) in which case ( |E.3| ) is replaced by expectation values in free theory and thus we 
find that 

a{e) = N^h{X,e/N^) (E.5) 

for some function h. 

As another example, let us look at the thermal average of ( |6.1C1| ), 

i J2 e-^^'^-e, = ^JdE e-^^niE) e{E) ^ e(i?^) (E.6) 

i 

Using the last equality of ( |6.10| ), the left hand side of ( [E.6| ) can in turn be written as 

Ep - {V)p (E.7) 

where (V) p is the thermal one-point function of V in the interacting theory and scales 
with N d.sO{N'^). Thus we can write 

e{E) = N'^g{\E/N'^) (E.8) 

for some function h. An exactly parallel argument yields 

E{e)=N^g{Ke/N^) . (E.9) 

To summarize, we find that the averaged values of r(e) and A(i?) are both of order 
0(N) in the 't Hooft limit for any nonzero A. Thus in the large limit, both the correlation 
length between interacting theory energy levels and the energy range that the free theory 
states are mixed under perturbation go to infinity. 
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Appendix F. Some useful relations 

In this appendix we derive some important relations which will be used in Appendix G 
to derive the matrix elements of an operator O between generic states in the high energy 
sector. 

F.l. Density of states 

The conservation of states implies that the density of states Q[E) of the full theory 
and Oo(e) of the free theory should be related by 

n{E) = Qoie{E))^§^ (F.l) 



which implies 



1 dn{E) _ de{E) 1 rfOo 

de{E) 



n{E) dE dE flo{e{E)) de 



In the large limit the second term in the above equation should be of order 0{1/N'^). 
Thus we find that 

I3{E) = Po{e{E))^§^ (F.2) 



with 



1 dn{E) 1 dQo .p.. 



We also expect that 

e(E{e)) ^ e (FA) 
Note that all the above relations are valid only to leading order in N. 

F.2. Properties of xe{^) and p^{E) 

Consider the microcanonical average of (|6.5| ) and which we denote as Pe{E) and 

Xe{^) respectively. From ( |6.5D and ( |6.6D one should have 



From (16.71) we should also have 

''dexE{e) = l (F.6) 
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and 

1 



dEp,{E) = ^ I dEn{E)xE{e) = 1 . (F.7) 

Given that 

e{E)= [ deexEie), E(i?) = A^(i?) = / de {e -e{E)yxE{e) (Fi 



we can write Xe{^) as 

XE{e) = fE{e-e{E)) (F.9) 

with fE a function which has a spread of A{E) ~ 0{N). Since we expect fEiyj) to faU off 
quickly to zero in the large limit outside the range (-iA(E), iA(S)), equations W^ ) 
and (^3) lead to 

/>oo 

dufE{u:) = l (F.IO) 



and 

0(E) 



dE-^fE{e-e{E)) = l (F.ll) 



Changing the integration variable of (|F.11| ) to e' = e(E) and using ( [F.4|) , we find that 



where in the second expression we have replaced /^(e') ■^^(e)" This is because, as a 
function of e — e', the spread of / is of order 0{N), while E{e') ^ E{e) + 0{^j^) ^ E{e). 
The second expression of ( [F.12|) can now be written as 

/oo 
dufE{uj)e-^°^^^^^^^ = 1 (F.13) 
-oo 

Equations (FUO|) and ( [FU3D can be written in a more symmetric manner as 



/oo rOO 
duje^^^^^'^ QEiuj) = / diue-^f^^^'^'^ gEiuj) = 1 (F.14) 
-oo J — oo 



where we have introduced a function 



Equations ([F.14|) imply that gE{i^) should fall off faster than e 2/^(^)1^1 as ct; ^ ±oo. 
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F. 3. A relation between matrix elements and correlation functions in free theory 

In this subsection we derive in free theory a relation between the matrix elements of 
an operator O between states in the high energy sector and correlation functions. For 
simplicity we consider theories with a single fundamental frequency ujq, like A/" = 4 SYM 
or (O). 



The Lehmann spectral decomposition for frequency space Wightman function G^^\uj) 
of some operator O in free theory can be written as 

(^) = ^ E e"^^Va6 S{u; -e, + e„) (F.16) 

a,b 

where pab = \ (fi|C'(0)|6) p. Due to energy conservation, O can only connect levels whose 
energy differences lie between —Aujq and Acuq, where A is the dimension of O, i.e. pat can 
only be non- vanishing for |ea — eb| < Awq- We can thus rewrite ( [I'MGI ) as 



1 A 

Gf{u) = ^ V Gkd{iv - kuo) (F.17) 

^0 



fc=-A 

with 



^0 V 



e-^'^Pkia) 



where ^^^^e denotes that one sums over \b) whose energy is given by = e. Note here 
we have assumed that the free theory energy levels are equally spaced as in A/" = 4 SYM 
theory on S^. We also introduced 

Pk{a)= E P-f' (^-19) 



We now separate the sum a in ( [F.18|) in terms of energies and degeneracies, i.e. 

E = EE 

We thus find that 

^^ = ^E^(^)^"'''^^(^) (F.20) 
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where we have introduced the micro-canonical average of pfc(6)for energy e 

We expect that the micro canonical average Pfc(e) should be a slow varying function of 
e, i.e. it can be written in a form N'^f{e/N'^) for some constant a. In the large 
limit since A/'(e)e~^'^ is sharply peaked at e/3 specified by, one can perform a saddle point 
approximation in (|F.20| ) to get 

Gk=M^p) + --- (F.22) 

^From ( [F.16| ) we thus find that 

G+iu) = Y^M^M^ - kuo) (F.23) 
k 

In the large limit since the connected part of G_(_(a;) scales with as 0{N^), thus we 
find from (g^) that 

p,(e^) ~ 0{N') . (F.24) 



Appendix G. Derivation of matrix elements 

In this appendix we give a derivation of ( |6.14|) . The main object of interests to us is 



El |2| |2 
|Cia| I'^jfel Pab 

a,b 

= El^*°I^E E 

a k ei,=ea+kiL!o 

Due to the sparse and random nature of pab, one cannot naively approximate the sums 
over a and b in by integrals. Instead one must be careful with the discreteness nature of 
the sum. Note that 

Y \Cjb\^Pab ^Cj{ea + k) ^ Pab = Cj{ea+k)pk{a) 

which can be justified as follows. Inside a given energy shell, pab can be treated as a 
random sparse matrix. Thus one can treat the summand as a random sampling of jcj^p. 
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Since the number of sampling points goes to infinity (as a power in A^) in the large 
limit, we can approximate |c,7,P by its average value of the energy shell. We now have 



Pij = X^X^ \cia\'^Cj{ea + k)pk{a) 

k a 

= ^^Cj{e + k)^\cia\'^Pk{a) (^q^-^ 



k e 



$^5^c,(e + /c)Ar(e)c,(e)p^(e) 



k e 

In the second line above we separated the sum over all states a into the sum over the 
energy and the sum over states in each energy shell. In the third line we replaced the sum 
in an energy shell by its average values. The replacement is all right since |ciap and pfc(a) 
are completely independent variables, so the average of their product should factorize. 

Now given that all quantities in the last line of ( |G.1| ) are averaged quantities, we 
approximate the sum over e by an integral. Averaging i,j over states of the same energy 
and using (|6.6D , we find 



/dt 
^^^^j^j^-^ XE2 (e + k)xE^ (e) Pk (e) 



5Z / + k) + ^' ~ - gi) Pfc(g) (G.2) 



k 



Oo(ei2+p + /c) 2 2 



with 



ei,2 = e(^i,2), €12 = \{e{Ei) + Z{E2)) = e{E), A12 = e{Ei) - e{E2' "^^^^ 



where 



E\ -\- E2 

E = , u! = E-i — Eo 

2 i ^ 



E 



Equation (|G.2|) can be further simplified as 

PEiE2 = Qj^Y'^^-^i^) J dP9E{p + k' + -uj)gE{p- - 



k 

1 



- -ujj 

(G.3) 

A{u;E) 



n{E) 
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with 

Gi2{k) = e-^^°(^-)'=p,(ei2), k' = -^k (G.4) 

dE 

and QEiyj) was defined in ( |i^'.15| ). Note that from equation ( |i^'.23D , G12 are essentially 
the Fourier components of free theory correlation functions. A[uj; E) should be a smooth 
function of u since the integral in ( |G.3|) appears to be well defined for all u. It is easy to 
check that 

A{-uj-E)=A{uj-E) (G.5) 
since Gi2{k) = Gi2{—k). Further as a; — » 00, we find that 

oce-i^(^)l"l . (G.6) 

Now let us examine possible singularities of A(a;, E) on the real axis. We start with 
the definition ( [F.9|) of /e- Since X-E(e) is the average of ( |6.6| ) over states of similar energies, 
fEiyj) must be a real positive function oi u E R. Then the function gsi^) introduced in 
( [b'.lSD should also be real and positive as e{E) is a monotonous function of E. The posi- 
tivity and normalization conditions (|F.14| ) imply that gsi'^) can at most have integrable 
singularities of the formEl 

9e{uj) ^ uj^uJi, ai<l (G.7) 

\UJ — UJi'°'i 



I 



Note that the closer ai is to one the smaller is Ki from the normalization requirement. 
Now let us look at the definition ( |G.3| ) of A{uj; E), 



A{uj- E) = J2 Gi2{k)s{uj + k) 
k 

where 

/■°° 1 1 

s{u) = I dxgE{x + -u)gE{x- -u) . (G.9) 



Note that the finite sum over k in (|G.8| ) cannot introduce singularities in u therefore we 
focus on s{uj). As gE{i^) falls off faster than e~5^^^^l'^l as u; — >• ±00, the integral in ( |G.9| ) 



Such integrable singularities can only arise if da have accumulation points in the large A'^ 
limit. While it appears unlikely that this can happen, we do not have a rigorous proof at the 
moment. 



48 



is convergent for a: ±00. Thus we only need to worry about possible divergences arising 
from the middle of the integration range. Integrating ( |G.9| ) we find that 

/oo 1*00 pOC 

dujs{uj) = / dxgE{x) / dygE{y) (G.IO) 
-00 J —00 J —00 

which is finite by ([F.14| ) . Therefore the only singularities allowed for s{ijj) are of inte- 
grable kind |^_^ |„ with a < 1. We can find the locations of ujs in terms of (integrable) 
singularities of gEi^^) as follows. Since gsix — ^u) and gEix + ^u) are both integrable the 
only possible divergences of (|G.9[) are at values of u for which the integrable singularities 



of two function sit on top of each other. This happens for u = Ui — ujj where the coi are 
the locations of the singularities for gE{oj). For uj = uji — ujj + e with e small the integral 
( |G.9| ) near x ^ + Wj) can be written as KiKj J^g dy |^_ 1 1 ^|c^ where S is some 

multiple of e. By rescaling we see that it behaves as . Therefore the integral s{uj) 

can at most have a singularity of the form |^_^'^^ |„ with a = ai + aj — 1 < 1. 

Thus we conclude that on the real axis A{u); E) can have at most integrable singular- 
ities of the form 

A{uj; E) oc as<l . (G.ll) 
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